Comment on "Critical temperature of trapped hard-sphere Bose gases" 
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In this comment, I discuss a recent path-integral Monte Carlo calculation by Pearson, Pang, 
and Chen ( Phys. Rev. A 58, 4796 (1998)). For bosons with a small hard-core interaction in a 
harmonic trap, the authors find a critical temperature which does not change with respect to the 
non-interacting gas. The calculation suffers from a serious discretization error of the many-particle 
density matrix. 
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In a recent article, Pearson, Pang, and Chen |ffl] studied 
hard-core bosons in a harmonic trap using Path-Integral 
Quantum Monte Carlo (PIMC) methods. The authors 
conclude that the critical temperature of this system for 
small hard-core diameters a does not notably change with 
respect to the non-interacting gas. These findings are dif- 
ficult to reconcile both with previous Monte Carlo work 
H and with mean-field calculations 0] , which agree very 
nicely with each other. However, the calculation of Pear- 
son, Pang, and Chen is flawed by a discretization error 
which affects the many-particle density matrix for small 
a. 

In 0], the density matrix p(R, R', (3) at inverse temper- 
ature (3 is reduced to an integral over high-temperature 
matrices p(R,R',t) with r = (3/M, and M the number 
of 'slices' 



(3) 



p(R,R',0) 



dR 2 



dR M p(R, #2,r) • ■ -p(R M ,R',T) 

(1) 



(cf Oj for notation, R = (ri, . . . , rjv) is the 3N- 
dimensional vector of particle coordinates, with N the 
number of particles; I use natural units with mass m = 1, 
oscillator length a — 1, h — 1). 

To evaluate p(Ri, Ri+i,r), the authors use the so- 
called "primitive approximation" (c/Q|) p(Ri, i?j+i,r) ~ 
pp(Ri,Ri+i,r) with 



Pp{R 1 ,R 1+ i,t) = cexp(- 



-Rj+i) 2 



2r 



-T^(i2i)) (2) 



For the present discussion, the harmonic potential is 
unessential and I will only consider the hard-core term. 

The interaction V(r) sets a small length scale, a, and 
the primitive approximation is justified only if, between 
Ri and Ri+i, the paths fluctuate by much less than this 
scale. Under this condition we can conclude that, if paths 
overlap neither at slice i nor at slice i + 1, they will not 
collide in between. The spatial fluctuations are deter- 
mined by the kinetic energy term in Eq. (^|), and one 
arrives at the condition 



In their work, Pearson, Pang and Chen consider, e. g. 
a value a = 0.02. It follows from Eq. (|J) that r must 
satisfy r < 0.0004. For a typical value of f3 = 0.2, this 
implies M > 500. In ref. Q, M - 10 is used. The ef- 
fect of the interaction is therefore largely underestimated 
and important deviations from the case of noninteracting 
bosons are not picked up. 

The condition r <C a 2 renders the use of the primitive 
approximation unpractical for large systems of weakly in- 
teracting bosons. However, it is very instructive to check 
the estimate Eq. (||) for a single pair of particles in the 
setting of ref. |Q. To do this, I have generated a large 
number of 3— dimensional Gaussian paths directly from 
the Levy construction (cf Q) for various values of M. 
These paths sample Eq. (0) for V = 0. p(R, R' , (3) is ob- 
tained by means of a correction factor, which differs from 
unity by the probability P of pairs of paths to collide at 
some slice i. 
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FIG. 1. Collision probability P for a pair of particles 
with hard-core diameter a = 0.02 as a function of their 
separation |r| (in reduced units, Oo = 1, /3 = 0.2), for 
M = 10, 100, . . . , 100, 000 (full lines, from below). Also shown 
is P in the limit M — > oo obtained from the eigenfunctions of 
the hard-sphere potential (dotted line). 
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As an example, I consider the diagonal density matrix 
for two particles separated by a distance \r\ (R = (0,r), 
R 1 = (0, r)). In fig. 1, the numerically determined proba- 
bility P is plotted for f3 = 0.2 as a function of |r| for values 
of r = 0.02 (M = 10) to r = 2 x 10~ 6 (M = 100,000). 
As predicted in Eq. (^), we need large values of M to 
obtain a converged result. In the limit M — > oo, P can 
also be computed directly from the eigenfunctions of the 
hard-sphere potential (cf Q), as presented in fig. 1. This 
shows that the interaction of two particles need not be 
calculated with Monte Carlo methods. This computation 
is at the heart of the recent successful PIMC calculations 
for weakly interacting bosons, both in the trap (|] || and 
in free space 0. 
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